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Abstract 

We examine the 1+1 dimensional CP{N) model in the large ./V limit by using 
the Schrodinger representation. Starting from the Hamiltonian analysis of the 
model, we present the variational gap equation resulting from the Gaussian 
trial wave functional. The renormalization of the theory is performed with 
insertion of mass and energy counter-terms, and the dynamical generation of 
mass and the energy eigenvalue are derived. 
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I. INTRODUCTION 



The Schrddinger representation [|IJJ2j has been proven to be a powerful method for prob- 
ing the non-perturbative aspects of quantum field theory It provides a natural way 
of extracting physical quantities as the eigenvalues of the equivalent operators acting on 
wave functionals, which describe the theory. An effective way of computing the Hamilto- 
nian eigenvalue is to use the variational method [[F|-f|. Especially, the Gaussian trial wave 
functional has been extensively used in a variety of models, and it reproduces the one-loop 



results 10 



In this paper, we reexamine the 1+1 dimensional CP(N) model |]TT]| for large N 



The purpose is to test the Gaussian approximation in a geometric theory with non-polynomial 
interactions, which is known to be exactly soluble in the large N limit. The study also gives 
some insight into gauge theories where similar features such as dynamical mass generation 
and asymptotic freedom arise. We approximate the vacuum wave functional with a Gaussian 
and calculate the energy eigenvalue by minimizing the expectation value of the Hamiltonian. 
The ensuing gap equation can be solved with a mass counter-term, and the solution exhibits 
the asymptotic freedom as expected. The renormalization of the Hamiltonian is performed 
with the insertion of an extra energy counter-term and the result is in agreement with the 
path integral treatment fT3|j . 

II. LARGE N CP(N) MODEL 



We first use the coadjoint orbit method to formulate the CP(N) model [14| in the large 
N limit. Let us introduce 

Q = gKg~\ (1) 

where g E G = SU(N + 1) and K = ^diag(A, -1, -1, .., -1). The action of the CP(N) 
model is given by 

S=-^-Jd 2 xTr(V tl QV»Q), (2) 

while A is coupling constant, and V M = d/dx^. Parameterizing g by (N + 1) vectors as 
g = (Zi, Z 2 , .., Zff+i) with Z p G C N+1 (p = 1,..,N + 1), such that Z* ■ Z q = 5 pq , and 
det(Zi, Z 2 , .., Z N+ i) = 1, we obtain 

S = I J d 2 x(V,Z{ ■ WZ, + {Z{ ■ V^) 2 ), (3) 

with the constraint Z\ ■ Z\ — 1, for Zx — (zi, ■ ■ ■ , z N+ i). This constraint can be solved by 
assuming z N+ i = z* N+1 . Introducing the Fubini-Study coordinate: 

^ = ViV^ , (4) 

ZN+l 

we obtain 



1 



Zn+i = i • (5) 



JN+\M 2 xlN + 

Substitution into (|3|) gives 

S = i J d 2 xg ab d^ a d^ b , (6) 
where g ab is the standard Fubini-Study metric on CP(N): 

( N + \^)5 ab - r*iJ b 2 2 2 

^ fe = (iV + 1^2)2 > M =\H +--' + Wn\ , (7) 

where tp a is an unconstrained iV-component field. Note that g ab can be obtained from the 
Kahler potential K(ip, = \n(N + |^| 2 ) by 

9ab d^dft ' 1 ' 

To perform Hamiltonian analysis, we introduce the conjugate momenta, 7r a and 7f a of ip a 
and tp a respectively by 

7T a = jg a bi> b , Tfa = ^fea^ 6 • (9) 

Then the Hamiltonian density 7i is given by 

H = nj a + nj- -£ = \g ab 7T a n b + I^V^V^ , (10) 
where g ab is the inverse of g ab : 

g ab = ^(N+\4j\ 2 )(N5 ab + r^ b ) . (11) 
The canonical quantization gives the following equal time commutators: 

[r(x),7r b (y)] = ih5 a b 5(x-y) , 

[r(x)Mv)] = ® a (x)My)] = o> (12) 

where x denotes only a space variable, and time is fixed at a common value. Using the 
differential representation of the momentum implied by the Schrodinger representation 

5 5 
the quantum Hamiltonian is given by (g = NX) 



2 



H 
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N 2 



dx 



-Ti\N + \M)(N5 ab + ip a ip 
N 3 (N+\^\ 2 )6 ab -^ a ij b 



2\2 



+ 



(14) 



As the CP(N) model generates a divergent mass in the quantum level, (see equation (|19D 
and (p0|) ), it is necessary to regularize by introducing a mass counter-term. In addition, the 
eigenvalue of the Hamiltonian turns out to diverge. Hence, we add the following part to the 
Hamiltonian 



H. 



c.t. 



dxji (N + 



+ £, 



C.t. 1 



(15) 



where and S c ,t. will be determined appropriately in the following in order to subtract the 
divergences of the Hamiltonian. Note that the mass counter-term of ( |I5D breaks the original 
SU(N + 1) invariance of CP(N) model, (|6]), to its linearly realized subgroup SU(N) x U(l). 



III. VARIATIONAL METHOD 



The Schrodinger equation will be studied within the variational approach [[7]J|. Let us 
take a Gaussian trial wave functional 



\l/[?/>, ijj] = exp < — 



dxdy(ip a (x) - ip a {x)) 



G ab (x,y) 
2h 



{My) -My)) 



(16) 



for some specific configuration ip, tf), and some propagator, G ah (x,y), to be determined. We 
suppressed the normalization constant. Then, the following expectation values result 



(4>a(x)) = 1p a (x) , (lp a { x )My)) = $a( x )4>b(v) + ^G ba (y, x) 



(ip a (x)ip b (y)) = Mx)My) > (M x )My)) = i a ( x )ib(y) > M x )n(y)) = jG b „(y,x) . (17) 

It is not possible to compute the expectation value of the Hamiltonian, (H), in a closed 
form, because of the non-polynomial character of H. Instead, we take the large N limit 
(with g fixed), and keep the dominant term in N. In order to simplify the result even more, 
we take ip to be x-independent. Due to the A r -plicity of the fields in the CP(N) model we 
can set the scales such that \ip\ 2 ~ N and G aa (x, x) ~ N. Hence, we obtain for large N 



(H) = | / d, 



- ^(N+ H 2 )(5(0)) 2 + \{N+ \j,\*)G£(x,x) + ^G aa (x,x)G^(x lX ) 



N 2 h 



g 2 (N+\ij\ 2 ) 

N 2 h 2 
' g 2 (N +\i)\ 2 ) 2 



V 2 x G aa (x,x') 



+ 



N 2 h 



2fc2 



g 2 (N + 

,2 N 



\2\2 



G aa (x,x) V 2 x G bb (x,x') 



(V x G aa (x,x')\ x , =x ) + -fr 2 (N+\i>\ 2 + hG aa (x,x))] +£ct. , (18) 
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where terms are kept only to order h 2 . In the following we are going to keep if) as an 
arbitrary parameter, while G^ 1 will be chosen such that it minimizes the expectation value 
of the Hamiltonian. The variation with respect to G a b(x,y) gives the following equation 

V N+\?p\ 2 J abK ' 

4 N 2 / hG cc (x,x)\/ „, 7 \ s . . . 

— s (1 v I J )( -Vl+m 2 )5(x -y)5 ab , 19 

where the finite mass, m 2 , is given by 



, = ("mih _ "G-^^-./^ 1( } s) s 2 s (20) 

TV V jv + |W;|2/ V4AT aaV 7 r ^° v 7 



This equation can be solved by performing Fourier transform of the delta function, 
resulting in 

G-J(x,y) = 2N s (l - ^££^1) f^J^^ e ^-y) 5ab (21) 



which yields the mass gap equation: 

V 2 = yJ ^P 2 + m2 ' ( 22 ) 

Note that the dynamically generated mass /x 2 has a quadratic divergence (times bare cou- 
pling) in terms of the momentum cutoff. This is consistent with the one-loop result in 
perturbation theory []. Substitution of the variational equation (|T9| ) into flUf) leads to effec- 
tive Hamiltonian. We can bring it into the following expression; 

In the above, we only kept terms up to order h for producing the one- loop result. In 
this expression the mass, m 2 , is finite. However, there are infinities connected with the 
ultra-violet behavior of the momentum. In the next section we are going to derive the 
renormalization of the physical quantities of the theory. 



IV. RENORMALIZATION 

The renormalization of the theory requires that the infinities appearing in the Hamilto- 
nian, are absorbed in the bare coupling, g, the mass, p,% and the energy counter-term, £ c t _. 
Introducing a momentum cut-off, A, the finite mass, m 2 , becomes 



1 Small ?/> a , expansion of the Lagrangian given in (||) contains the derivative interaction of the 
form —l/^^dip^ldipl 2 + \$d , ip\ 2 ). This lowest order interaction will produce quadratic divergence 
in the one- loop order contributing to the dynamically generated mass. 



4 



m 



2 , 2 
/i + /i 



t4 



9 / A 2 . m2 . 2i 2As 

— (A + — + m log — 

47r v 2 m ' 



(24) 



where the divergences for A — > oo, have to be absorbed in //q and g. We have set h — 1. In 
order to treat the infinities, we first define a finite renormalized coupling constant as 

1 1 1 , M 

with an arbitrary renormalization scale M. We also define a renormalized mass by 



In terms of (|2oT) and 



^ l-^] o- 

]), relation ( p^ ) becomes 



(26) 



m 



2 , 9R m 



(l + log 



M 2 



(27) 



The definitions ( p5|) and (|26|) can also be viewed as determining //q and (7 for finite values of 
[i 2 R and with respect to the cut-off, A. Eq. (|25D shows asymptotic freedom as expected. 
In terms of the renormalized coupling the Hamiltonian in (|23|) becomes 



w = J v l y (fa| (_ + _, og _ + 



m 
8^ 



(28) 



where we have chosen the counter-term £ c t such that it cancels the appearing quadratic 
divergence, j-A 2 , of the Hamiltonian. This renormalized Hamiltonian is not unique because 



the final expression is dependent upon the renormalization prescription in (25) . If one 
chooses another definition by adding a constant C/An term to the right-hand side of (p5|) , 
the final Hamiltonian fl2|) changes by A{H) = N[f dx](-Cm 2 /8ir). Note that relation (|27|) 
remains the same. The constant C can be fixed by requiring a renormalization condition. 
If we demand, for example, 



d(H) 



dm 2 



where ml is the solution of 



M 2 exp[ 



0. 



8vr 
9r 



(29) 



(30) 



we find C = 2, and this (H) agrees with the one in the literature |13]| . Note that (H)\ m 2 =m 2 
does not depend on the coupling constant explicitly, a phenomenon known as dimensional 
transmutation 11131. 
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V. CONCLUSIONS 



We have shown that the variational technique with Gaussian wave functional reproduces 
the known results of CP(N) model in the large N limit. Unlike the path integral method 
which uses an auxiliary field, Hamiltonian procedure necessitates two counter terms. One 
is a mass counter-term to treat the divergence appearing in the mass gap equation and 
to extract a finite mass. The other is a quadratically divergent energy counter-term. The 
remaining divergences in the Hamiltonian can be absorbed by renormalizing the coupling 
constant and mass. 

The present work can be extended in a couple of ways. First, recall that nonlinear sigma 
model with target manifold of symmetric space can also be solved in the large N limit 
375 1 . It would be worthwhile to test the Gaussian method in this case also. Secondly, it 



would be interesting to perform higher-loop analysis which requires a systematic expansion 
of Gaussian wave functional JTB]. 
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